Average observational quantities in the timescape cosmology 
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We examine the properties of a recently proposed observationally viable alternative to homo- 
geneous cosmology with smooth dark energy, the timescape cosmology. In the timescape model 
cosmic acceleration is realized as an apparent effect related to the calibration of clocks and rods of 
observers in bound systems relative to volume-average observers in an inhomogeneous geometry in 
ordinary general relativity. The model is based on an exact solution to a Buchert average of the 
Einstein equations with backreaction. The present paper examines a number of observational tests 
which will enable the timescape model to be distinguished from homogeneous cosmologies with a 
cosmological constant or other smooth dark energy, in current and future generations of dark energy 
experiments. Predictions are presented for: comoving distance measures; H[z); the equivalent of 
the dark energy equation of state, w{z); the Om{z) measure of Sahni, Shafieloo and Starobinsky; 
the Alcock-Paczyhski test; the baryon acoustic oscillation measure, Dy-, the inhomogeneity test of 
Clarkson, Bassett and Lu; and the time drift of cosmological redshifts. Where possible, the predic- 
tions are compared to recent independent studies of similar measures in homogeneous cosmologies 
with dark energy. Three separate tests with indications of results in possible tension with the ACDM 
model are found to be consistent with the expectations of the timescape cosmology. 

PACS numbers: 98.80.-k 98.80.Es 95.36. +x 98.80.Jk 



I. INTRODUCTION 



The paradigm for our current standard model of the 
universe assumes that the universe is well-described by 
a geometry which is exactly homogeneous and isotropic, 
with additional Newtonian perturbations. The under- 
lying geometry is assumed to be that of a Friedmann- 
Lemaitre-Robertson- Walker (FLRW) geometry, and in 
matching the cosmological observables that derive from 
such a geometry, we have been led to the conclusion over 
the past decade that the present-day universe is domi- 
nated by a cosmological constant or other fluid-like "dark 
energy" with an equation of state, P — wp, which vio- 
lates the strong energy condition. 

Although the matter distribution was certainly very 
homogeneous at the epoch of last-scattering when the 
cosmic microwave background (CMB) radiation was laid 
down, however, in the intervening aeons the matter dis- 
tribution has become very inhomogeneous through the 
growth of structure. Large scale surveys reveal the 
present epoch universe to possess a cosmic web of struc- 
ture, dominated in volume by voids, with galaxy clusters 
strung in sheets and filaments that surround the voids, 
and thread them. Statistical homogeneity of this struc- 
ture appears only to be reached by averaging on scales 
of order 100/i~^ Mpc or more, where h is the dimen- 
sionless parameter related to the Hubble constant by 
= 100/ikm sec~^ Mpc""'^. The problem of fitting a 
smooth geometry to a universe with such a lumpy mat- 



ter distribution [H, Q is a nontrivial one, but central to 
relating observations to the numerical values of the av- 
eraged parameters which describe the Universe and its 
evolution as a whole. 

Given the observed inhomogeneity of the present epoch 
universe, a number of cosmologists have questioned 
whether the FLRW geometries are adequate as a descrip- 
tion of the universe at late times 0-0- In particular, 
the deduction that the universe is accelerating might in 
fact be a result of trying to fit the wrong cosmological 
model. One central question in the fitting problem is the 
issue of deriving the average evolution of the inhomoge- 
neous geometry. If one considers irrotational dust cos- 
mologies, and averages just inhomogeneous scalar quan- 
tities, in Buchert's scheme [Sj one finds an average of 
the Einstein equation in which there is a Friedmann-like 
evolution modified by backreaction [l^, [H, 113 , 
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where an overdot denotes a time derivative for ob- 
servers comoving with the dust of density p, a{t) = 

[V(t)/V(to)]^^^ with V{t) = d^x^/det^, angle brack- 
ets denote the spatial volume average of a quantity, so 

that (7^) = (^J^d^x^/det^n{t,x.)^ /V{t) is the average 

spatial curvature, and 

- ((O- W')-2(a2), (5) 
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is the kinematic backreaction, = ^aapcy"^^ being the 
scalar shear. We use units in which c = 1. Eq. ^ is 
an integrabihty condition needed to ensure that eq. ([1]) 
is the integral of eq. 

One must be careful in interpreting equations ([T])~ 
([5]) since the spatial averages refer to average quantities 
which depend on the domain of integration on a spatial 
hypersurface. Observers measure invariants of the local 
metric, not a spatially averaged metric, and cosmological 
information comes to us on null geodesies. Given these 
problems, the Buchert approach has been criticised [13], 
and the whole area of backreaction is the subject of some 
debate and controversy. In recent work [2^ [2^, HO] I have 
developed an interpretation of solutions to the Buchert 
equations which circumvents the criticisms of Sec. 3 of 
ref. [l^l • It differs from other approaches to the Buchert 
equations that have been used in the literature [3H-[36j. 
As well as circumventing objections that have been raised 
against Buchert averaging, the new interpretation has a 
conceptual basis which can be understood as an extension 
of the equivalence principle , and it leads to a quanti- 
tative model universe with predictions [28, 38] which thus 
far are in good agreement with observation. In particu- 
lar, by Bayesian comparison the RiessO? gold supernovae 
la (Snela) data set [39| agrees with the model predictions 
at a level which is statistically indistinguishable from the 
standard spatially flat ACDM [H, [40|] . The same best-fit 
parameters also fit the angular scale of the sound hori- 
zon seen in CMB data, and the effective comoving baryon 
acoustic oscillation (BAO) scale seen in ang ular diameter 
tests of galaxy clustering statistics [H, 

Given these promising indications, it is important that 
the cosmology of refs. j28l [29j is developed well beyond 
the stage of what might be regarded as a "toy model" , so 
that it can be confronted by all the same observational 
tests that are applied to the ACDM model. For example, 
many current precision tests involve the detailed analysis 
of the CMB [4J|, and of galaxy clustering statistics [45j- 
[s^ l. To construct tests of similar precision will require 
the development of new numerical codes for the analy- 
sis of large datasets adapted to the present cosmology, 
analogous to those based on the decades of detailed work 
that have been applied to the standard cosmology. 

Such goals represent an arduous project, and here I 
will simply take a few steps in the direction of confronting 
the observations. The aim of the present paper is not to 
present a detailed analysis of current data sets, but to 
outline a number of observable quantities which might 
be tested in future. Since the predictions obtained for a 
number of these quantities can be readily compared to 
existing independent analyses of homogeneous cosmolo- 
gies with dark energy, I will make relevant comparisons 
where possible. I will confine the discussion here to aver- 
age quantities which are relevant at all redshifts on scales 
greater than the scale of statistical homogeneity. Other 
relatively local tests which deal with quantities within 
the 100h~^ Mpc scale of statistical homogeneity (28l. Issj 
will be left to future work. 



The plan of the papers is as follows. In Sec. |TI] I will 
summarize the key features of the model introduced in 
refs. [H, [2§|, while also providing some further discus- 
sion. Additional technical details which were not pro- 
vided in ref. [29.] on account of space restrictions are given 
in the Appendices. In Sec. lIIIl I discuss the luminosity dis- 
tance and angular diameter distance relations, and their 
interpretation in terms of the equivalent of a "dark energy 
equation of state" , which enables a direct comparison to 
recent studies to be made. In Sec.|TV]related diagnostics, 
H{z) and the Om{z) measure are evaluated and discussed 
in relation to recent studies. The Alcock-Paczyhski and 
BAO tests are treated similarly in Sec. |Vl The expected 
nontrivial signature of a test of the Friedmann equation 
of Clarkson, Bassett and Lu is determined in Sec. IVII A 
prediction for the Sandage-Loeb test of the time drift of 
cosmological redshifts is presented in Sec. IVIIl Sec. IVIIII 
contains a concluding discussion. 

II. OVERVIEW OF THE TIMESCAPE MODEL 
A. Voids and walls 

I will begin by briefly reviewing the two-scale "frac- 
tal bubble model" |28| which I am hereby renaming the 
"timescape model" [5J] , concentrating on the operational 
interpretation of observations. The model is constructed 
by identifying the observed scales most relevant to the 
observed present epoch inhomogeneous structure as be- 
ing negatively curved voids and spatially flat walls, which 
surround bound structures. 

Since galaxies and galaxy clusters formed from per- 
turbations which were greater than critical density, then 
given an observable universe which on average has nega- 
tive Ricci scalar curvature, we have a natural separation 
between walls and voids. As there is assumed to be a gra- 
dient in spatial curvature, it is assumed we can always 
enclose the bound structures which formed from over- 
critical perturbations within regions which are spatially 
flat on average, and marginally expanding at the bound- 
ary. These boundaries are called finite infinity regions 
[28|, with local average metric 

ds^ = -dr^ + aJ{T^) [dijl + rj^dn^] . (6) 

The walls constitute the union of such finite infinity re- 
gions. Observationally they would correspond to all ex- 
tended structures that contain galaxy clusters, namely 
sheets, filaments and knots [55*1 . 

Voids of a characteristic diameter 3Qh^^ Mpc are ob- 
served to fill 40-50% of the universe at the present epoch 
[5^. In addition there are numerous minivoids, which 
have been well-studied in the local volume [53]. To- 
gether voids of all sizes appear to dominate the volume 
of the present epoch universe, the exact volume fraction 
depending on the empirical definition of a void in terms 
of some particular negative density contrast. In the two- 
scale approximation of refs. [1^ [2^ both the dominant 
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voids and minivoids are assumed to be characterized by 
the same negatively spatial curvature scale, with a metric 
at the void centres being given by 

= + fd^'- + smh'irj,)dfl^] (7) 

We construct an average over the disjoint union of wall 
and void regions over the entire present horizon volume 
V = VjCi'^, where 

= /viOv^ + /wiflw^, (8) 

/vi and /wi = 1 — /vi being the respective initial void 
and wall volume fractions at last scattering. The finite 
infinity scale only becomes operationally defined once 
regions start collapsing and structure forms. Thus at 
last scattering a different interpretation of the wall and 
void components in ([H]) is required. At this epoch the 
wall fraction is understood as that fraction of the 
present horizon volume which comprises perturbations 
whose combined mean density is the same as the mean 
density of the statistical ensemble of perturbations, in- 
cluding those beyond the horizon. The void fraction, /v, 
is understood to be that (small) fraction of the present 
horizon volume in underdense perturbations which was 
not compensated by overdense perturbations at last scat- 
tering. It is convenient to rewrite ^ as fv{t) + fw{t) = 1, 
where fw{t) = fwiO-w /o,^ is the wall volume fraction and 
fvit) = fvidv /a^ is the void volume fraction. 

B. The scale of "statistical homogeneity" 

In order to physically identify observables we must 
specify what is to be identified as a "particle" of dust. 
The question of what constitutes a particle of dust is not 
directly addressed in Buchert's scheme, although perhaps 
implicitly many researchers think of galaxies as being the 
particles of dust, as historically this is the way the matter 
is treated in the FLRW model. However, galaxies evolve 
considerably over time and are not homogeneously dis- 
tributed at the present epoch. Thus if we wish to follow 
cosmic evolution from the epoch of last scattering to the 
present, with no assumptions about homogeneity, then 
we must coarse grain the dust on scales over which mass 
fiows can be neglected, so that each dust particle remains 
of a roughly fixed mass, even if the mass differs somewhat 
from particle to particle. 

Here we take a "dust particle" to be of at least the 
scale of statistical homogeneity, 100h~'^ Mpc or some- 
what larger |_58j . The scale of statistical homogeneity is 
taken to refer to a scale volume within which the struc- 
ture of voids and walls is roughly similar, if such a box 
is chosen at random on a spatial slice in the observable 
universe at late epochs. It is important to realise such 
volumes will not have the same density. Rather they will 
have a density which is distributed about a mean with a 
standard deviation of order several percent, by an argu- 
ment that follows from eq. ^ below. 



We must stress that the universe is not considered to 
be a FLRW model, and our terminology of a "scale of 
statistical homogeneity" is not the same as in a FLRW 
model. A scale of homogeneity, in its sense in the FLRW 
model, is not assumed to exist. The principal difference 
is that in a FLRW model the density of the observable 
universe, sampled over the present horizon volume, is 
assumed to be the mean density of the ensemble from 
which our observable horizon volume was drawn. In the 
FLRW case the standard deviation of the density of spa- 
tial volumes would decrease to zero as one sampled ever 
larger volumes greater than the homogeneity scale, as is 
the case for any stationary stochastic process (59| . 

Such a state of affairs cannot be expected to prevail, 
however, given cosmic variance and an initial spectrum 
of density contrasts of all possible length scales which are 
nested within each other, which is the expectation from 
primordial inflation. Given cosmic variance, then as one 
samples larger and larger volumes that become compa- 
rable with the horizon volume, one is dealing with fewer 
and fewer individual fluctuations rather than a statisti- 
cal ensemble. The assumption in the timescape scenario 
is that the present horizon volume is underdense relative 
to the ensemble mean density, which at last-scattering is 
extremely close to critical. 

The fact that it does nonetheless make sense to think 
of a "scale of statistical homogeneity" as above, however, 
is a simple consequence of the fact that although the 
density perturbations have all possible length scales, the 
magnitude of these contrasts was strongly bounded at the 
time of last scattering. In other words, given a universe 
which was close to homogeneous at last scattering, it can 
only evolve so far from homogeneity within the finite age 
of the universe. 

The relevant scale for a cutoff to the "scale of statisti- 
cal homogeneity" is scale of the largest acoustic wave in 
the plasma at last scattering - of order llQh^^ Mpc. The 
simple reason for such a cutoff, is that below this scale ini- 
tial density contrasts may be amplified by acoustic waves 
in the plasma, so that rather than having initial density 
contrasts of say 5p/p ^10~^ in nonbaryonic dark mat- 
ter, the initial density contrasts will be somewhat larger. 
The second acoustic peak in the CMB anisotropy spec- 
trum - i.e., the first refraction peak - for example, will 
amplify the density contrast of underdense regions, and 
may therefore be the feature of the primordial spectrum 
responsible for the fact that the dominant void fraction 
is associated to a specific scale 30h^^ Mpc [56j . 

Given that some initial density perturbations are am- 
plified below the acoustic scale, and that the CMB 
anisotropy spectrum is fairly flat at long wavelengths, 
the acoustic scale provides a cutoff analogous the cut- 
off between the nonlinear and linear regimes of struc- 
ture formation, although here there is no single global 
FLRW model about which such regimes linear are de- 
fined [iOl- Below the scale of statistical homogeneity we 
will typically find density contrasts ~ 1 which char- 

acterize the nonlinear regime, as is the observed case for 
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30h~^ Mpc diameter voids [s^. Above the acoustic scale, 
we can be sure that the perturbations at last scattering 
have very similar amplitudes as a function of scale. Al- 
though the perturbations in the photon-baryon plasma 
have contrasts 6p/p^ 10^^ at this epoch, the density con- 
trast in nonbaryonic dark matter is expected to be some- 
what larger; e.g., of order S p / p ^ 1Q~'^-10~^ , depending 
on one's dark matter model. 

The standard deviation of the density of cells on scales 
larger than the scale of statistical homogeneity can be es- 
timated crudely by assuming that such cells each evolve 
as an independent Friedmann universe from a smooth 
perturbation at the epoch of last scattering. This approx- 
imation is justified since the relevant scale is the one over 
which there are no appreciable average mass flows from 
one dust cell to another. We assume that the backreac- 
tion contributions do not dominate the volume-average 
evolution, and make our rough estimate from the Fried- 
mann equation with prcssurelcss dust only, for which 
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This leads to a present epoch density contrast 



Spt 



(9) 



where the density contrast is relative to the critical den- 
sity, so that dpt — flj^jit) — 1 etc, where f2^^ is a den- 
sity parameter for the isolated region only. (Physically, 
the critical density is that within a spatially flat wall re- 
gion.) Thus if we take Spt ~ 10^'* at last scattering, when 
z ~ 1090 and when H ~ 2/{3t) with t ~ 380, 000 yr, we 
are led to dp^ ~ 0.025/h^ ~ 0.06 if /i ~ 0.65. 

This crude estimate can be compared to the actual 
density variance determined from large scale structure 
surveys [ill, HI]. Sylos Labini et al [6l| have recently 
determined the variance in the number density of lumi- 
nous red galaxies (LRGs) in the SDSS-DR7 by dividing 
the full sample of 53,066 galaxies in the redshift range 
lO^'* < z < 0.3 into N equal nonoverlapping volumes. 
Over the range 4 < < 15, the standard deviation is 
found to be of order 8%, consistent with an earlier mea- 
surement of 7% by Hogg et al [6l[ in a smaller LRG 
sample. These values are very close to our order of mag- 
nitude estimate of 6%. Provided LRGs are correlated 
to the actual density, then the variance in the percentage 
density contrast will be commensurate. In fact, such vari- 
ances can be used to constrain the dark matter density 
contrast at last scattering. A measurement of 8% would 
indicate, reversing the argument above, that a contrast 
of 5p/p^\Q^^ in nonbaryonic dark matter at last scat- 
tering is an order of magnitude too large. 

Given a nearly scale-invariant spectrum of density per- 
turbations, with perturbations nested in perturbations, 
our expectation is that the variance in density should 
not decrease appreciably if sample volumes are increased 
at nearby redshifts. In principle, it should be possible to 
calculate it as a function of scale, given the constraints 



from the CMB anisotropy spectrum at long wavelengths. 
For spatial slices at higher redshifts, looking further back 
in time, the variance would decrease in accord with ^ - 
provided that a sample of objects such as LRGs can be 
found which does not exhibit strong evolutionary effects 
over the range of redshifts in question. 

Of course, the estimate based on ^ could be further 
refined to take backreaction into account; but further 
accuracy can only be gained when one has a tighter esti- 
mate of the dark matter density contrast than simply an 
order of magnitude. Furthermore, the statistical physics 
of cosmic structures in the timescape scenario may well 
differ from that of the FLRW model [5^ in significant 
ways; one has to revisit the whole problem from first 
principles. 

In summary, the observed universe is not assumed to 
be homogeneous or to approach any single global FLRW 
model at any scale. The "statistical homogeneity scale" 
- which will coincide roughly with the BAO scale - repre- 
sents a scale above which the variance in density contrasts 
is bounded at the 10% level, and below which density 
constrasts become as large as they can possibly be. 



C. The bare Hubble flow and bare cosmological 
parameters 

Given our identification of dust particles coarse- 
grained at the scale of statistical homogeneity, and pos- 
sible very large differences in spatial curvature and grav- 
itational energy within such a cell, we do not assume 
that Buchert average time parameter, i, is the relevant 
parameter measured by every isotropic observer - those 
who see an isotropic CMB - within any dust cell. Rather 
it is the time parameter measured by an isotropic ob- 
server whose local spatial curvature happens to coincide 
with the Buchert volume-average spatial curvature (7^) . 
We employ an ansatz of an underlying quasilocal uniform 
Hubble flow within a dust cell, below the scale of statis- 
tical homogeneity, in terms of local proper lengths with 
respect to local proper times, which both vary with gradi- 
ents in spatial curvature and gravitational energy. This 
ansatz provides an implicit resolution of the Sandage- 
de Vaucouleurs paradox [1^, and can be understood in 
terms of a generalization of the equivalence principle [STj . 

The metrics ^ and ([7]) are assumed to represent the 
local geometry for isotropic observers at finite infinity 
and at void centres respectively. Within the scale of 
statistical homogeneity the metrics ([6]) and ([7]) are as- 
sumed to be patched together with a condition of uniform 
quasilocal bare Hubble flow [H, [s^l 



H 



drw 



dflv 

dTv 



(10) 



which will preserve isotropy of the CMB. The mean CMB 
temperature and angular anisotropy scale will vary with 
the gradients in gravitational energy and spatial curva- 
ture, however. 
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For the purpose of the Buchert average we refer all 
quantities to one set of volume-average clocks: those that 
keep the time parameter t of eqs. H])-© so that 



a 



where 



1 da^ 



and 



1 dov 
flv dt 



(11) 



(12) 



and 



7,, 



dt 
dr™ 



and 



7v 



dt 

dTv ' 



(13) 



are lapse functions of volume-average time, t, relative to 
wall and void-centre observers respectively. The ratio of 
the relative Hubble rates hr = H^/Hy < 1 is related to 
the wall lapse function by 



7, 



= 1 



(1 - hr)f.. 



(14) 



and 7^ = ^r7w 

The Buchert equations for pressureless dust with 
volume-average density p^^ are solved p9l | in the two- 
scale approximation by assuming that there is no back- 
reaction within walls and voids separately [b^, but only 
in the combined average. With this assumption, the kine- 
matic backreaction term becomes ,28J 



Q - 6/v(i - /v) (Hv - H^r - 



2/v 



3/v(l - /v) 



(15) 



The resulting independent Buchert equations consist of 
two coupled nonlinear ordinary differential equations [1^ 
for a{t) and fv(t), which may be written as 



M 



n, 



= 0, 



(16) 
(17) 



Equations ([T6l) , p?|) are readily integrated to yield an 
exact solution |29(| , which is listed in Appendix [X] to- 
gether with its simple tracking limit in Appendix [b1 For 
initial conditions at last scattering consistent with obser- 
vations, solutions are found to reach within 1% of the 
tracking limit by a redshift z ~ 37 29] . Thus the tracker 
solution will be used for the purposes of the specific cos- 
mological tests which are investigated in this paper [65| . 

It should be noted that for the solution found in ref. 
[2^, the backreaction term is at most of order 4.2% [66l |. 
Its redshift dependence for the best-fit parameters is ex- 
hibited in Fig. [1] Although ^Iq is negative, it is never 

large enough relative to fj^^ to dominate the r.h.s. of 
(PT|) and give volume average cosmic acceleration. The 
backreaction itself is not the sole reason for apparent cos- 
mic acceleration; that is also a question of how volume- 
average evolution is interpreted in terms of a local metric. 




FIG. 1: The bare backreaction density parameter CIq as a 
function of redshift for the timescape model with /vo = 0.762, 
= 61.7 im sec'^ Mpc"^ . 



where 
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2/3 f 1/3 



9/v(l - /v)-ff' 



(18) 
(19) 
(20) 



are the volume-average or "bare" matter density, curva- 
ture density and kinematic backreaction density param- 
eters respectively, and p^^^ being the present epoch 
values of a and p^. The average curvature is due to 
the voids only, which are assumed to have fcv < 0. The 
volume-average deceleration parameter is given by 



2"m 



(21) 



D. Dressed cosmological parameters 

One must take care in physically interpreting the solu- 
tion of the Buchert equation, since it does not represent a 
single exact solution of Einstein's equations, but rather a 
spatial average. Observers measure invariants of the local 
metric, and information carried by radial null geodesies 
from distant parts of the universe. In ref. ^281 a means 
of interpreting the Buchert equation was developed as 
follows. 

First, since cosmological information is obtained by 
a radial spherically symmetric average, we construct a 
spherically symmetric geometry relative to an observer 
who measures volume-average time, and with a spatial 
volume scaling as a'^(i). 



ds^ 



-df + a^{t) df + A{n,t) dif 



(22) 
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where the area quantity, A{'r],t), 



satisfies 



a'(0Vi(^„)/(4^), fy„ being the 



conformal distance to the particle horizon relative to 
an observer at 77 = 0, since we have chosen the particle 
horizon as the scale of averaging. The metric is 
spherically symmetric by construction, but is not a 
Lemaitre-Tolman-Bondi (LTB) solution since it is not 
an exact solution of Einstein's equations, but rather of 
the Buchert average of the Einstein equations. 

In terms of the wall time, r„, of finite infinity observers 
the metric is 



-7w(T-w)dr2 + d^{T^) df + A{f], Tw) dn^ . (23) 



However, this geometry, which has negative spatial cur- 
vature is not the locally measured geometry at finite in- 
finity, which is given instead by Since ^ is not a 
global geometry, we match (jH]) to ([^5]) to obtain a dressed 
wall geometry, which is effectively the closest thing there 
is to a FLRW geometry adapted to the rods and clocks 
of wall observers. The matching is achieved in two steps. 
First we conformally match radial null geodesies of (O 
and (j23p . bearing in mind that null geodesies are unaf- 
fected by an overall conformal scaling. This leads to a 
relation 



dr]^ 



7w(l-/v)'/' 



(24) 



along the geodesies. Second, we account for volume and 
area factors by taking 7/^ in ([6]) to be given by the integral 

of HMD. 

The wall geometry ([S]), which may also be written 



.2/3 



[dr;2 + jf^dn" 



(25) 



on account of ([5]), is a local geometry only valid in spa- 
tially flat wall regions. We now use pi]) and its integral 
to extend this metric beyond the wall regions to obtain 
the dressed global metric 



ds' = -d7 



■df 



-a" (1 - /v) 



2/3 



^2 ' .2/3 

/ W 7 W" 



■?7w('7,T-w)df^^ 
= ~dTl + a'{T^)[df + rl{n,T^)dn'] (26) 
where a = 7^^flj and 

Whereas ([6]) represents a local geometry only valid in spa- 
tially flat wall regions, the dressed geometry ([26|l extends 
as an average effective geometry [67,] to the cosmologi- 
cal scales parameterized by the volume-average confor- 
mal time, which satisfies d^ = dt/a = dr^/a. Since the 
geometry on cosmological scales does not have constant 
Gaussian curvature the average metric (pS)) . like (|22p . is 
spherically symmetric but not homogeneous. 



In trying to fit a FLRW model to the universe, the cos- 
mological parameters we obtain effectively have numeri- 
cal values close to those of the dressed geometry (|26l) . In 
particular, we infer a dressed matter density parameter 



M 



a dressed Hubble parameter 



H = 



1 da 
a dTw 



1 da 
a dTw 



7w 



(27) 



(28) 



and similarly a dressed deceleration parameter, where the 
overdot still denotes a derivate with respect to volume- 
average time. As demonstrated in refs. [H, [2^ in a void- 
dominated universe the dressed deceleration parameter is 
negative at late epochs, even though the bare decelera- 
tion parameter (j2ip is positive. Thus cosmic acceleration 
is realised as an apparent effect due to the variance of lo- 
cal geometry from the average, leading to variance in the 
calibration of clocks and rods. 

In the rest of the paper we will drop the subscript "w" 
from both Tw and 7^, as we will not need to make explicit 
reference to the time measured in void centres. Thus t 
and 7 will be assumed to refer to wall time. 



III. COMOVING DISTANCE D{z) AND 
EQUIVALENT OF THE "EQUATION OF STATE" 

In testing fluid-like dark energy scenarios, or modified 
gravity theories that can be cast as an effective fluid with 
an equation of state P = wp, a common question is: how 
can the equation of state parameter, w{z), be constrained 
as a function of redshift? Unfortunately, if dark energy 
is some purely unknown physics, then it is completely 
unclear how one should expand it as a power series. A 
linear series in z will not converge for z > 1, for example, 
so series in z/(l 4- z) are sometimes considered. Unless 
one has a precise physical model of dark energy to be 
tested, then any constraints are completely dependent 
on how one chooses to characterize such a power series. 
When constraints on the value of w from cosmological 
observations are quoted in the literature, it is often on 
the basis that w is simply a constant, even though there is 
no known physics for making such an assumption, apart 
from the cosmological constant case of w = — 1. 

In this section I will derive the equivalent of the "equa- 
tion of state" style observational tests, although the ter- 
minology "equation of state" does not have a meaning in 
terms of actual observables, given that the model in ques- 
tion is not characterized by a fluid with = wpj^. Let 
us recall that in the case of the standard FLRW models, 
the equation of continuity for such a dark energy com- 
ponent in a background universe with scale factor a(i), 
viz.. 



Pd 



3^(1 
a 



)Pz,=0, 



(29) 
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may be integrated to give 

using a^/a — 1 + z, where it is assumed that the equa- 
tion of state parameter varies with redshift. To obtain 
an expression for the luminosity distance one substitutes 



(I30p in the spatially flat Friedmann equation for matter 
plus dark energy, 




and uses the resulting expression for d, to determine the 
conformal time integral 



' FLRW 



dt 



% da 



f7,,„(l + z')3 + f^^„exp 3/„^ 



(l+tti(z"))dz" 
1+z" 



(32) 



1 



AW 



where ^^^^q 

nosity distance is then given by c?^ 
The quantity 



D 



and 



DO 



The standard lumi- 
,(1 + ^). 



1 + z 



(33) 



is the comoving distance quantity directly related to the 
luminosity distance. The angular diameter distance is 
also related by 



-4 1 + z 



{l + zf 



(34) 



We observe from (|32p that H^D docs not depend on the 
value of the Hubble constant, H^, but only directly on 
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Given observed quantities such as the apparent 
luminosity-redshift relation or an angular size-redshift 
relation for standard candles or standard rulers, we can 
take derivatives of (l32l) to obtain 



w{z) 



{l + z)D'-^D" + 1 



z)^H'^D'^ - 1 



(35) 



where the prime denotes a derivative with respect to z. 
This gives a formal equation of state to any comoving dis- 
tance relation, assuming an underlying spatially flat dark 
energy model. Such a relation can be applied to observed 
distance measurements, regardless of whether the under- 
lying cosmology has dark energy or not. We should note, 
however, that such a w{z) has first and second derivatives 
of the observed quantities, and so is much more difficult 
to determine observationally than direct fits to a quantity 
such as D{z). 

For the timescape universe, equivalent comoving, an- 
gular diameter and luminosity distances can be defined 
in terms of the dressed geometry ([6]). We have a dressed 
luminosity distance relation 



ap(l + z)r^ 



(36) 



where a^^ — 7g ^Oq, and the effective comoving distance 
to a redshift z is D ^ a^r^, where 



7-w 7 (1 - /v) 



1/3 



d<' 



7(t')(l-/v(i'))^/^«(i') 



As discussed in Sec. [TTl since spatial sections are not of 
constant Gaussian curvature this effective comoving dis- 
tance represents a fit to our spatially flat rods once radial 
null geodesies are conformally matched, and geometric 
factors are taken into account. 

For the tracker solution (|Bip . (|B2p the cosmological 
redshift satisfies 



z + 1 



where 



an7 



(2 + /v)/v'/' _ 2W/3(t + 6) 



a% Vli'^H^t /;/^5„i(2i + 35)4/3 



1/3, 



^_ 2(l-/vo)(2 + /vo) 
The integral in ([57]) is readily evaluated to give 



(38) 
(39) 



^ 1 + z 



2dt' 



t {2 + Mt')){tr/' 

t^/^{Ht,)-m) (40) 



where 



- 2t + g ln(^^2/3_;,l/3il/3 + 52/3 
61/3 /2tl/3_5l/3 

H ^ tan ;= — ■ 

VS V \/36i/3 



(41) 



It is straightforward now to compare distance mea- 
surements in the timescape model with those in spa- 
tially flat ACDM models. The timescape model which 
best fits the Riess07 gold data set had a void fraction 
/vo = 0.76lQ;Qg, and dressed Hubble constant = 
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61.7^]^'^ km sec ^ Mpc ^, where la uncertainties are 
quoted In Fig. [U we plot 



//„i? = if//3[^(ij-^(t)](l + z) 



(42) 



for the best-fit model with /vo = 0.762, as compared to 
three spatially flat ACDM models with different values 
of ^j^jQ, (or of il^g = 1 — ^j^jq)- Fig. O shows that over 
redshifts between the present epoch and last scattering, 
the timescape model interpolates between ACDM mod- 
els with different values of ^j^^q- For redshifts z < 1.5 
Urpg is very close to ^^qdj^j for the parameter values 
(f}^^Q,17^o) = (0.34,0.66) (model (iii)) which best-fit the 
RiessOT Snela data only [Sgl- For very large redshifts 
that approach the surface of last scattering, z < 1100, 
on the other hand, D^^ very closely matches -D^^^j^ 
for the parameter values (fi^j^g, fi^p) = (0.249, 0.751) 
(model (i)) which best-fit WMAP5 only [ij]. Over red- 
shifts 2 < z < 10, at which scales independent tests are 
conceivable, D^^^ makes a transition over corresponding 
curves of D^^-,j^j^ with intermediate values of (f^^^Q, ^^q)- 



The D 



ACDM 



curve for joint best fit parameters to Snela, 



BAO measurements and WMAP5 ;44l, ($1^^^^, fJ^g) = 
(0.279,0.721) is best-matched over the range 5 < z < 6, 
for example. 

Given the difference of D^^ from any single ^'^^dj^ 
curve becomes pronounced only in the range 2 < z < 6, 
it may be difficult to distinguish the models on the basis 
of the measurement of alone from BAO surveys, which 



will be able to measure d^(z) up to 1% to z < 3. How- 
ever, joint measurements of other parameters, such as 
H{z), may make for definitive tests, as will be discussed 
later. Gamma-ray bursters (GRBs) do probe distances 
to redshifts z < 8.3, and could be very useful. There has 
already been much work deriving Hubble diagrams using 
GRBs. (See, e.g., [1^.) It would appear that more work 
needs to be done to nail down systematic uncertainties, 
but GRBs may provide a definitive test in future. An 
analysis of the timescape model Hubble diagram using 
69 GRBs has just been performed by Schaefer [t^I, who 
finds that it fits the data better than the concordance 
ACDM model, but not yet by a huge margin. As more 
data is accumulated, it should become possible to distin- 
guish the models. 



A. Effective "dark energy equation of state" 

The equivalent of an equation of state, w{z), for the 
timescape model may be determined from and ([M)) . 
The specific analytic expressions for the first and second 
derivatives of D are 



dD 
d^ 



t 



{2t -b)d, + {2t + 36)' 



3 (2i2 + m + 252) 



(43) 
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3.5 
3 

2.5 
2 

1.5 
1 

0.5 







(b) 



(i) 

(ii) 

(iii) 



200 400 



600 
z 



800 1000 



FIG. 2: The effective comoving distance H^D{z) is plotted for the best-fit timescape model, with /vo ~ 0.762, (solid line); and 
for various spatially flat ACDM models (dashed red lines). The parameters for the dashed lines are (i) = 0.249 (best-fit 
to WMAP5 only); (ii) fi^^^ = 0.279 (joint best-fit to Snela, BAO and WMAP5); (iii) 0.^^^^ = 0.34 (best-fit to RiessO? Snela 
only). Panel (a) shows the redshift range z < 6, with an inset for z < 1.5, which is the range tested by current Snela data. 
Panel (b) shows the range z < 1100 up to the surface of last scattering, tested by WMAP5. 
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d^D _-t {2t + 36) [2 it + b) {2t + 56) {2t^ + 3bt - 6^) + (2t + 36) (St" + 266^^ + 53bH^ + 56bH + 186^)] 



dz2 



9(1 + z) (2i2 + 36i + 262)-' 



(44) 



SSI) and dMl) in ([35]) and use the fact that by jBgJ, = (4/2^ + /vo + 4)i7g/[2(2 + /vo)], to obtain 



In these expressions d. is given by (|40p and t is given impHcitly in terms of the redshift, z, via (I38p . We now substitute 



(45) 



(40t^ - 286^4 - 27462i3 _ 34963i2 _ Q2bH + 246^) d^ +t{2t + 3bf {20t^ + 56bt^ + mH - 46^ 
{(2t - 6) + (2t + 36)^} |^o(z + 1)3*2 (2t - 6)d^ + (2t + 36)^ ^ - 9 (2*2 + 36i + 262) 



where 







»Afo(4/v^O + /vO + 4)^ 

4(2 + ur 



(46) 



and is given by ((39|) and ([40)1 . In fact, fi^^^ 
i(l-/vo)(2 + /vo), so that 

(l-/vo)(4/v2o + /vO + 4)2 



(a) 



8(2 + /vo) 



-1.5- 




(b) 



FIG. 3: The artiGcial equivalent of an equation of state i)45l) . 
constructed using the effective comoving distance i)35|). ptot- 
ted for the timescape tracker solution with best-fit value 
fvO = 0.762, and two different values of ^j^.^y (a) the canon- 
ical dressed value fi^^-, 
"mo = 0.279. 



i(l-/vo)(2 + /vo) = 0.33; (b) 



Since the w{z) expression is an artificial mathematical 
construction for the present model, we can also determine 
w{z) if a value of i^^^Q different from the canonical value 
i(l — /vo)(2 + /vo) is assumed. In this way, we arrive 
at the example w{z) curves plotted in Fig. [3] The fact 
that the denominator of (|45p goes through zero means 
that w{z) becomes formally infinite and changes sign at 
a value of z which depends on the value of f^^^g assumed 
fn\ . This feature illustrates how pointless it is to talk 
about an equation of state of dark energy, or to choose 
to "reconstruct" w{z) if the underlying unknown physics 
has nothing to do with a fluid in the vacuum of space. 
What is actually measured is a quantity such as D{z), 
illustrated in Fig. [21 and this is perfectly smooth. 

Phenomenologically, for the canonical best-fit dressed 
value of = 0.33 [3^, one finds that w{Q) ~ -0.758 
and that w{z) crosses the "phantom divide" w{z) ~ —1 
at z ~ 0.464. The average value of u;(z) ~ —1 on the 
range z < 0.7, while the average value of w{z) < —1 if 
the range of redshifts is extended to higher values. This 
agrees with the evidence of the Snela data. 

In fact, in a recent study [TSj which examines con- 
straints on the equation of state by combining the Con- 
stitution Snela data 'i^ with WMAP5 ;44l and SDSS 
constraints, Zhao and Zhang find 95% confidence level ev- 
idence in favour of a model with w{z) > — 1 e (0.25, 0.5), 
w{z) < — 1 G (0.5,0.75), meaning that w{z) must cross 
the phantom divide in the range 0.25 < w < 0.75. The 
fiducial model of Fig. ^a) crosses the phantom divide 
almost in the centre of this redshift range. 

Another recent investigation [t^I draws different con- 
clusions about evidence for dynamical dark energy. How- 
ever, while the results of Serra et al [73| are consistent 
with a cosmological constant at the 2a level, they are also 
consistent with the best-fit timescape model at the same 
level, as illustrated in Fig. [H The different conclusions 
drawn by the authors of refs. ,73] and [74j| result not only 
from using somewhat different data sets, but also from 
differences in the treatments of data bins. 

In considering Fig. [J] it should also be borne in mind 
that there are significant systematic issues between the 
SALT and MLCS data reduction methods, as will be dis- 
cussed in Sec. IVIIIl The Union [4^ and Constitution 
compilations use the SALT method. A new analysis 
of 103 SDSS-II Snela in the redshift range, 0.04 < 
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B. Angular— size redshift relation 




FIG. 4: The artiEcial equivalent of an equation of state i|45l) 
is compared with a recent analysis of Serra et al Jt^] . The 
third panel of Fig. 1 of ref. J74j is combined with the curve 
of w{z) for the best-Rt value /vo ~ 0.76;too9 (solid curve; la 
limits dotted curves). Following ref. [74] 2a uncertainties are 
plotted. The la uncertainties are tabulated in Table\^ 



redshift 








0.0 
0.25 
0.5 
0.75 
1.0 


-0.97 ±0.22 
-1.05 ±0.10 

-0 71+0-44 

1 79+0.73 
^- ' ^-0.81 


-0.86 ±0.13 
-1.04 ±0.07 

-1 06+°''^ 

-0.471^3^ 


-1.02 ±0.28 

1 O1+0.49 
— i.»»_2.76 



TABLE L Values of w{z) determined by Serra et al [74j us- 
ing a standard FLRW cosmology are compared to the artifi- 
cial equivalent of w{z) for the timescape model: and w^-, 
are the values determined by combining WMAP5 CMB and 
SDSS-DR7 BAO data with the Union and Constitution Snela 
data sets respectively, as given in Table I of ref. [tJ]. The 
equivalent Wrpg for the timescape model uses a void fraction 
fvO = 0.761q Jg as determined in ref. ^3§\. la uncertainties 
are listed in each case. 



z < 0.42 when combined with 185 Snela from other sur- 
veys, yields best-fit parameters w — — 0.96±0.06 (syst)± 
0.12 (stat) and Qj^^^ = 0.265 ± 0.016 (syst) ± 0.025 (stat) 
using SALT-II to fit to spatiaUy flat FLRW models with 
constant w, but w = —0.76 ± 0.07 (syst) ± 0.11 (stat) and 
= 0.307±0.019 (syst)±0.023 (stat) using MLCS2k2. 
Use of MLCS data reduction is therefore also likely to 
somewhat change the data values in Table [J and Fig. [D 

At this stage the uncertainties, especially systematic 
ones in data reduction, are too large to draw firm conclu- 
sions, but future measurements may change the picture. 
Of course given specific models of dark energy, greater 
statistical leverage is obtained simply by comparing H^D 
directly on a model by model basis. 



The angular size, S — i/d^, of a class of objects of 
uniform proper length, is readily determined from (j40p . 
(HT|) . Empirically the differences from the ACDM model 
are not very large. For the best-fit value /vo = 0.762 the 
minimum angle occurs at z = 1.74, as opposed to z = 
1.67 for a spatially flat ACDM model with Q^^^ = 0.249, 
or z = 1.56 for a spatially flat ACDM model with ^j^jq — 
0.34 (see Fig.[5]). The angle subtended by standard rulers 
in the timescape model is very slightly less than for the 
comparison spatially flat ACDM models. At z = 6 the 
difference is of order 9-15%. 




FIG. 5: The angle, 5 (in arcsec), subtended by a lOkpc 
source as a function of redshift for the timescape model with 
fvO = 0.762, = 61.7 im sec~^ Mpc~^ (solid line) as com- 
pared to the equivalent angular size relation for three spa- 
tially Bat ACDM models (dashed lines from top to bottom) 
(a) (!r2^^Q,n^(,) = (0.279,0.721), = 71.9 im sec"^ Mpc"^ 



(b) (n 

(c) (n 



MO ' ^AO ) 
AfO' ^Ao) 



: (0.249, 0.751), = 71.9 im sec-^ Mpc- 
(0.34, 0.66), = 62.7im sec'^ Mpc'^ . 



IV. THE H(z) AND Om(z) MEASURES 

Recently Sahni, Shafieloo and Starobinsky [t^ have 
proposed a new diagnostic of dark energy [771 1 the func- 
tion 



Om{z) — 







- 1 



(1 



(47) 



on account of the fact that it is equal to the constant 
present epoch matter density parameter, f^^^Q, at all red- 
shifts for a spatially flat FLRW model with pressureless 
dust and a cosmological constant, but is not constant if 
the cosmological constant is replaced by other forms of 
dark energy. For a spatially flat universe with pressure- 
less dust plus some arbitrary dark energy one has 



Om(z) = r?,, ± (1 - n.. ) 



(1 + ^) 



3(l+u)) 



- 1 



(l + zf 



(48) 
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For general FLRW models H = D'-^ + Vl^^H'^D'^ 

only involves a single derivative of D{z), and so the diag- 
nostic (|47)) is easier to reconstruct observationally than 
the equation of state parameter, w{z). 




FIG. 6: The dark energy diagnostic Om{z) of Sahni, ShaBeloo 
and Starobinsky Jfdij plotted for the timescape tracker solu- 
tion with best-Et value /vo = 0.762 (solid line), and la lim- 
its (dashed lines) from ref. \3Hl: (a) for the redshift range 
Q < z < 1.6 as shown in ref. \8ul] : (b) for the redshift range 
0<z <6. 

The quantity Om(z) is readily calculated for the 
timescape model, and is shown in Fig.[6l What is striking 
about Fig. [SI as compared to the curves for quintessence 
and phantom dark energy models plotted in ref. [76|, is 
that the z = intercept 



Om(0) - f H'\, = 



3/vo 



■4)(2 + /vo) 



(4/v'o + /vO + 4)2 



(49) 



is substantially larger than in the dark energy mod- 
els. We note from p9|) that lim/^p_»o O"i(0) — 1, and 



lim* 



Oto(O) = |, with a minimum value of Om{0) ~ 
0.638 at =i 0.774. The best-fit present epoch void 
fraction [3^| gives a value of Oto(O) very_close to the 
minimum. For the range /vo = 0.76tQ g^ [38] Om(0) is 
tightly constrained to the range 0.638 < Om{Q) < 0.646. 



A further difference for the timescape model is that 
Om{z) does not asymptote to the dressed density param- 
eter ^j^jQ in any redshift range. For quintessence mod- 
els Om{z) > i^j^jQ, while for phantom models Om{z) < 
il^^p, and in both cases Om{z) — > f^^^Q as z ^ oo. In 
the timescape model, Om{z) > ^j^j^^ — 0.33 for z < 1.7, 
while Om{z) < ^j^jq for z > 1.7. It thus behaves more 
like a quintessence model for low z, in accordance with 
Fig. [3] However, the steeper slope and the completely 
different behaviour at large z mean the diagnostic is gen- 
erally very different to that of dark energy models. For 
large z, 



lim Om{z) ~ 



2(l-Ao)(2 + /vo)' 



(4/^0 + /vo 



4)2 



(50) 



For 



giving a value i^j^jQ < Om{co) < i^^jQ, if /vO > 0.25. 
example for /vo = 0.762, we find Om(oo) ~ 0.2. 

Shafieloo, Sahni and Starobinsky [80| have recently 
tested the Om{z) statistic against CMB, BAO and Snela 
data, including the Constitution Snela data 'i?] . In com- 
paring their results with Fig. [6] it should be noted that 
their analysis entails taking particular empirical func- 
tions for w(z), and then best-fitting the free parame- 
ters. The two functions they choose are: (i) w{z) = 
Wi^-\-w^z/{l-\-z); and (ii) w{z) — — i[l-|-tanh((z — zt)A)], 
where w^, Wj^, Zt and A are empirically fit constants. 
In both these cases w{z) is monotonic and cannot com- 
pletely accommodate the equivalent "artificial dark en- 
ergy equation of state" for the timescape model as de- 
picted in Fig. [3] at large z. Furthermore, the effec- 
tive w{z) of Fig. [31[a) becomes nonlinear in the range 
0-5 ^ z < 1, contradicting the parameterization of case 

(i) of 80]. Also, it crosses the "phantom divide" at 
z ~ 0.464 contradicting the parameterization of case 

(ii) . However, we can expect that the empirical forms 
of w{z) assumed by Shafieloo et al to have some compar- 
ative value for the timescape model at very low values 
of z. The greatest leverage should come as z ^ 0. It 
is therefore interesting to note that of the two empirical 
forms for w(z) assumed by Shafieloo et al, the one that 
provides the better fit, case (ii), gives a best-fit intercept 
Om(0) remarkably close to the expectation from (I49p for 
/^o = 0.762, viz. Oto(O) = 0.638. (See Fig. 3 right hand 
panel of ref. [13] ■) Since this is not the expectation for 
either a typical quintessence or phantom energy model, it 
is an encouraging result for the timescape model, which 
is also consistent with the study of Zhao and Zhang [t^ . 

Shafieloo et al have suggested [80] that their analysis 
of the recent data might give a hint that "dark energy is 
decaying" . Given that the results of ref. [s^ appear to be 
consistent with the expectations of the timescape model, 
our analysis sheds a different light on this interpretation. 
It should also be noted that a pure FLRW model with 
substantial negative spatial curvature, i.e., with ^^^^ > 0, 

whose value 



will give an intercept Om(0) = ilj^^^ + |^1j.q, 
could assume a similar value to that obtained for the 
timescape model. Of course, this would require a value 
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of rjj,p which is ruled out by the WMAP analysis for the 
FLRW case, which is why such values have not been con- 
sidered by Shafieloo et al. As observed above Om(0) has 
a very tight range of values for a wide range of reason- 
able values of /vo- Thus if the tests of the Om{z) statistic 
could be improved to include a wider range of empirical 
functions, including those that more closely mimic 
our relation then this would be an interesting test 
once significantly more data becomes available. 
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FIG. 7: The function H~^H{z) for the timescape model with 
fvO ~ 0.762 (solid line) is compared to H~^H{z) for three spa- 
tially flat ACDM models with the same values of {^j^jq, ^^ao) 
as in Fig.\^ (dashed lines): (i) J^^^) = (0.249,0.751); 

(ii) (fi^o'^Ao) = (0.279,0.721); (iii) 
(0.34,0.66). 



The strong differences seen in the Om{z) diagnostic be- 
tween the timescape model and typical dark energy mod- 
els might be seen to arise from the fact that it accentuates 
the differences which already exist in the dressed H(z) 
function, which is quite different to that of the Friedmann 
equation. Using (jB9|l we plot H{z)/Hg for the best-fit 
timescape model, and compare it to the spatially flat 
ACDM models that were plotted in Fig. H For z < 1.5 
H{z)/H^^ for the timescape model with f^o — 0.762 is 
greater than for the ACDM models shown. The abso- 
lute value of H{z) is partly compensated for, however, 
by the higher value of that is generally assumed for 
the ACDM models. 

Gaztafiaga, Cabre and Hui 81] have recently given 
measurements of H{z) at three redshifts, inferred from 
the separation of radial and transverse BAO scales in 
the SDSS-DR6 data, as will be discussed in Sec.|Vl How- 
ever, their values are model dependent, being estimated 
according to 



Hiz)t,, 



BAO 



i?ov/0.25(l + z)3-|-0.75 



' WMAP 



with a fiducial expansion rate for a spatially flat ACDM 
model, with i^^g = 0.25, used to convert redshifts to 
distances. Any estimates of H(z) will inevitably involve 
some model dependence, unless one can perform a test 



such as the time drift of cosmological redshifts, which 
will be discussed in Sec. IVIII 



V. THE ALCOCK-PACZYNSKI TEST AND 
BARYON ACOUSTIC OSCILLATIONS 

Alcock and Paczynski devised a test [sj] which relies on 
comparing the radial and transverse proper length scales 
of spherical standard volumes comoving with the Hubble 
flow [s^. This test was originally conceived to distin- 
guish FLRW models with a cosmological constant from 
those without a A term. The test is free from many evo- 
lutionary effects, but relies on one being able to remove 
systematic distortions due to peculiar velocities. 

For the timescape model the Alcock-Paczyhski test 
function determined from the dressed geometry is 



1 



'AP 



6z 



HD _'i {2t^ + 'M + 29) (1 + z)d^ 



t(2t- 



3bfz 



(51) 

where t is given implicitly in terms of z by (|38p. 

In Fig. [8] the Alcock-Paczynski test function (|5ip is 
compared to that of spatially flat ACDM model with dif- 
ferent values of (flj^j^,fl^^). The curve for the timescape 
model has a distinctly different shape to those of the 
ACDM models, being convex. However, the extent to 
which the curves can be reliably distinguished would re- 
quire detailed analysis based on the precision attainable 
with any particular experiment. 

Current detections of the BAO scale in clustering 
statistics of LRGs [i^-fEOl can in fact be viewed as a 
variant of the Alcock-Paczyhski test, as they make use 
of both the transverse and radial dilations of the fiducial 
comoving BAO scale to present a measure 



zD^ 



1/3 



Df 



-1/3 
AP ■ 



(52) 



In Fig. [9] the BAO radial test function ([52]) is compared 
to that the same spatially fiat ACDM models plotted in 
Fig. [51 for the same redshift ranges. 

Although the Dy measure for the timescape model is 
significantly different to that of the ACDM model at the 
higher redshifts shown in Fig. [9jb), we see from Fig.[9l[a) 
that at the nearby redshifts the Dy measure gives consid- 
erably less discriminatory leverage. A case in point is pro- 
vided by the ratio J'g 35.9 2 = (0.35)/!),^ (0.2), which 

In this case 

0.35:0.2 ~ 



has been determined observationally [46, 5 
the timescape model with /vo = 0.76lg'Qg gives 
1.632j:°:°°f, as compared to r^^^,^^ = 1.664+°;°°? for a 
spatially fiat ACDM model with f^^^^ = 0.28 ± 0.03; val- 
ues which are close. By comparison the observed ratio 
was initially estimated to be gg.g ^ = 1.812±0.060 [46j . 
but this estimate has recently been revised to g^.g ^ = 
1.736 ±0.065 [53. 

In fact, one must exercise caution in comparing the 
prediction of rg gg.g ^ for the timescape model with the 
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FIG. 8: The Alcock-Paczynski test function /^p = 7 1 1| | 
for the timescape model with /vo = 0.762 (sohd hue) is 
compared to f^p for three spatially Bat ACDM models 
with the same values of (Slj^^g, 57^^) as in Fig. [2| (dashed 
lines): = (0.249,0.751); (ii) (0^0- ^^ao) = 

(0.279,0.721); (iii) (J^^^^, fi^J = (0.34,0.66). Two redshift 
ranges are shown: (a) < 2: < 1; (b) < z < 6. 



"observed" ratio [H, [s^] since the galaxy clustering data 
has been analysed in a manner which assumes an un- 
derlying FLRW model. The relevant analyses [i^ [lOl in- 
volve transformations to Fourier space to treat the power 
spectrum. To revisit such an analysis for the timescape 
model is far from trivial, as it requires a recalibration of 
transfer functions, and of the cosmological drag epoch, 
Zd, for a rather different cosmological parameterization 
given that we are dealing with a model which does not 
evolve as a homogeneous isotropic cosmology. In par- 
ticular, the mass ratio of nonbaryonic dark matter to 
baryonic matter can be somewhat different [ssj from the 
concordance ACDM model, and this needs to be consid- 
ered. Given that the difference in the value of r„ „_ „ „ 

U.oo:U.z 

quoted between refs. 46] and [50] is due to changes in 
the manner in which the data is treated, as well as the 
fact that there is more data, it is clear the differences in 
calibration due to a change of the nonbaryonic to bary- 
onic mass ratio could also similarly affect the value of the 
"observed" ratio. 




(b^ 1 2 3 4 ^56 

FIG. 9; The BAO radial test function H^D^ = H^Df"^''^ 
for the timescape model with /vo ~ 0.762 (solid line) is 
compared to H^Dy for three spatially Hat ACDM models 
with the same values of (r2^^g,f2^g) as in Fig. [2| (dashed 
lines): (i) (n^,„,n^„) = (0.249,0.751); (ii) (^^^o'^ao) = 
(0.279,0.721); (iii) (J^^^^, = (0.34,0.66). Two redshift 

ranges are shown: (a) < 2 < 1; (b) < z < 6. 



A derivation of tools which would enable us to per- 
form the BAO tests to the extent of refs. [i^, HO] is be- 
yond the scope of the present paper. Instead, it is our 
aim to simply explore what the best possible discrimi- 
nating tests will be. In this regard, we note that if we 
compare Fig. [5fa) and Fig. [l^a), then it is clear that 
the Alcock-Paczyhski test provides much more signifi- 
cant differences between the timescape model and ACDM 
models than the Dy measure. In fact, the Dy measure 
is currently employed because there is not yet sufficient 
data to separately estimate both the radial and trans- 
verse BAO signals directly, as would be required for the 
Alcock-Paczyhski test. 

Gaztahaga, Cabre and Hui [8l| have recently claimed 
to separate the radial and angular scales corresponding to 
the BAO in the 2-point correlation function, by assum- 
ing a nonlinear gravitational lensing magnification bias. 
Using SDSS-DR6 data they have exhibited a correlation 
function in both the radial and transverse dimensions, 
for redshift shces at z = 0.15-0.30 and at z = 0.40-0.47. 
They have not yet provided separate estimates of both 
the radial and transverse BAO scales. However, pro- 
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redshift 
range 




"so 




0.15-0.30 


0.132 


0.028 


3.7 


0.15-0.47 


0.12 


0.026 


3.6 


0.40-0.47 


0.124 


0.04 


2.1 



TABLE II: Values of O^^^h^ 



'-'^ inferred by Gaztanaga, 
Cabre and Hui [81(, and the resulting mass ratio of nonbary- 
onic dark matter to baryonic matter, fl^Q/flg^. 



vided their techniques are robust, then a direct Alcock- 
Paczynski test may soon be on the horizon. Naturally 
such estimates will have model dependence. From the 
the point of view of the timescape model, one must care- 
fully consider not only the treatment of redshift space 
distortions, but also any assumptions which rely on cali- 
brations of FLRW models, as discussed above. 

One point of the analysis of Gaztanaga, Cabre and Hui 
is suggestive. They find some tension between their best- 
fit value of the baryon density parameter 



BO 



0.06 and 

the WMAP5 value [44| fl^^ ~ 0.0432. The discrepancy 
is greater in the higher redshift slice. Their results are 
summarized in Table [Til The inferred values for the mass 
ratio of nonbaryonic dark matter to baryonic matter of 
ripp/ri^P ~ (^A/o ~ ^so)/^BO "^"^ whole sam- 

ple, 3.7 in the lower z slice and 2.1 in the higher z slice, 
as compared to the expectation of a ratio of 6.1 from 



WMAP5, for which (r2„n/i 



(0.0227,0.1308) 

|44l |. In other words, the best-fit values indicate a some- 
what higher mass fraction of baryons than the fit to 
WMAP5 with a FLRW model. This is confirmed by 
analysis of the 3-point correlation function [sHi , and is a 
feature which the authors find difficult to explain as a sys- 
tematic error. The analysis of the 3-point function yields 
a best fit m n^^jg = 0.28 ± 0.05, l^^^, = 0.079 ± 0.025. 

For the timescape model by comparison, analysis of 
the Riess07 gold data [H, Ha] yields dressed parame- 
ters Qj^j^ = 0.33tg:ig, flg^ = 0.080;Ej;EJi3, and a ratio 
flu^/flgg = 3.II24 from supernovae alone. Demand- 
ing a fit of the angular diameter distance of the sound 
horizon [281 to within 4% would reduce these bounds 



to il^^/flgg = 3.1_j;3 for the timescape model. Thus 
the higher baryon density indicated by the analysis of 
Gaztafiaga, Cabre and Hui is consistent with the expec- 
tations of the timescape model. 

Finally, we note that although the reality of the BAO 
measure is accepted by most researchers [ill -[HI, Sylos 
Labini et al ^3] have questioned this. Although Sylos 
Labini et al detect the BAO scale in the LRG sample, 
they point out that its amplitude is less than the over- 
all density variations of 8% at large scales, and further- 
more the correlation function remains positive where the 
ACDM model predicts it should be negative. Sample un- 
certainties may limit the strength of this conclusion [5l[ , 
however. 

In my view, although the results of ref. [62] may poten- 



tially indicate problems with a statistical analysis based 
on the expectations of a FLRW cosmology, the BAO is 
a real feature which will survive despite the observed in- 
homogeneities. The point is that given a universe which 
was very close to homogeneous and isotropic at last scat- 
tering, it can only evolve so far away from homogeneity in 
the time available since that epoch. Thus there is every 
reason to expect that statistical analyses of the type that 
are being performed can pick up a feature in the two- 
point correlation function, even if there are larger scale 
variations in density of order 8%. The exact properties 
of the statistical correlation functions within a frame- 
work such as the timescape cosmology await a detailed 
analysis. The main difference is that the density of the 
observable universe when measured on scales larger than 
that of "statistical homogeneity" will retain some intrin- 
sic variance, and furthermore is not the time-evolution 
of the mean density of the statistical ensemble at last 
scattering. This is likely to have important consequences 
for the statistical analysis. 



VI. TEST OF (IN)HOMOGENEITY 

Recently Clarkson, Bassett and Lu [s^l have con- 
structed what they call a "test of the Copernican prin- 
ciple" based on the observation that for homogeneous, 
isotropic models which obey the Friedmann equation, the 
present epoch curvature parameter, a constant, may be 
written as 



kO 



[Hiz)D'{z)]^ 1 



(53) 



for all z, irrespective of the dark energy model or any 
other model parameters. Consequently, taking a further 
derivative, the quantity 

C{z) = l + H^{DD" -D'^) + HH'DD' (54) 

must be zero for all redshifts for any FLRW geometry. 

A deviation of C{z) from zero, or of (|53p from a con- 
stant value, would therefore mean that the assumption 
of homogeneity is violated. Clarkson, Bassett and Lu re- 
fer to this as a "violation of the Copernican principle" . 
Given the viewpoint outlined in ref. [28j . simply associ- 
ating FLRW models with the Copernican principle is too 
great a restriction on its general philosophy. One should 
distinguish the Copernican Principle, which is generally 
understood as the statement that we do not occupy a 
privileged position in the universe, from the Cosmologi- 
cal Principle that the universe is described by a spatially 
homogeneous isotropic geometry. 

In the presence of inhomogeneity there can still be sta- 
tistically average cells - taken here to be of size lOOh^^ 
Mpc - but with a variance of the geometry within such 
cells. As observers in an average galaxy, our position is 
unremarkable from the point of view of the Copernican 
principle. Nonetheless, the local geometry in an aver- 
age void can be markedly different from the geometry in 
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an average galaxy. Given that observers and the things 
they observe are necessarily in bound structures, struc- 
ture formation provides a selection effect in terms of our 
local geometry vis-a-vis the volume-average geometry in 
a void. Given this improved understanding of the Coper- 
nican principle, one should not call the test of Clarkson, 
Bassett and Lu a test of the Copernican principle. It is 
simply a test of the validity of the FLRW models. 

Since the timescape model is inhomogeneous, it will 
certainly violate the test of Clarkson, Bassett and Lu. 
If one can determine H{z) in a model independent way, 
then tests of relations ([55]) or (|54p could not only rule 
on whether the FLRW model is violated, but also test 
the timescape model. Analytic expressions for HD' and 
HD" are obtained by multiplying and dSJ) by ((B8|l . 
Combining the results with and (jBlip we find that 
(l53l) becomes 



B{z) 



(55) 



where 



B = 



{2t - b)d^ 



{2t + 3bY 



{2t - b)d^ 
(2< + 3bf 



while (|54p becomes 



{2t-b)d^ 36(10^2 + ll&t- 262)^2 
{2t + 36)^ 



t{2t + 'iby 



(57) 



Once 



again, is given by (j40p and t is given implicitly 
in terms of z via (|38p . 

We plot the functions B{z) and C{z) in Figs. [TUl and 
nn The function C differs appreciably from the FLRW 
value of zero. However, two derivatives are required to 
determine C{z), which is subject to greater uncertainties 
for actual data, given that D{z) is effectively what is 
measured. Thus it would be more feasible to determine 
B{z) — [HD'Y — 1, which involves a single derivative of 
the observed curve. (It makes more sense to plot B{z), 
rather than the right hand side of eq. (j53p . which in- 
volves a division by zero as 2; ^ 0.) In Fig. [Tn]S(z) for 
the timescape model is compared to the expectations for 
ACDM models with a small amount of spatial curvature, 
as compatible with WMAP. The form of B{z) is very dif- 
ferent at small redshifts, which suggests that this will be 
a useful observational test. Furthermore, since B{z) has 
a maximum value and also changes sign, its form for the 
timescape model is very different to that of any FLRW 
model. In the FLRW case B{z) is always a monotonic 
bmction whose sign is determined by that of ^^^.q. At 
large z, or equivalently at early times ast^O,;B(z)— s-0 
and C{z) — *■ for the timescape model, consistent with 
the fact that it coincides with a spatially flat Einstein-de 
Sitter universe at early times. Since C{z) involves second 
derivatives, it goes to zero more slowly than B{z): for the 
best-fit solution of Figs. [TO] and [IH S(llOO) ~ -0.0029, 
while C(llOO) ~ 0.075. 



It is interesting to compare Fig. [TT] with the corre- 
sponding plot of C{z) for a LTB model with a large void 
recently given in Fig. 14 of ref. The magnitude of 

C{z) is considerably larger for the timescape model. 




FIG. 10: The homogeneity test function B{z) = [HD']'^ - 
1 is plotted for the timescape tracker solution with best-fit 
value /vo = 0.762 (solid line), and compared to the equivalent 
curves B = Q.^^^{H^Df for two different ACDM models with 



small curvature: 



(56) (b) n = 0.28, n 



(a) n^^^ = 0.28, n^, 
= 0.73 and Q^n 



0.71 andO., 



■ 0.01; 



-0.01. A spatially 



flat FLRW model would have B{z) 




-0.2- 



-0.3 ^ 



FIG. 11: The homogeneity test function C{z) given by i|54D is 
plotted for the timescape tracker solution with best-fit value 
/vo = 0.762. Any FLRW model would have C{z) = 0, regard- 
less of its spatial curvature. 



VII. TIME DRIFT OF COSMOLOGICAL 
REDSHIFTS 

For the purpose of the (in) homogeneity test consid- 
ered in the last section, H{z) must be observationally 
determined, and this is difficult to achieve in a model 
independent way. There is one way of achieving this, 
however, namely by measuring the time variation of the 
redshifts of different sources over a sufficiently long time 
interval 89], as has been discussed recently in relation 
to tests of (in)homogeneity by Uzan, Clarkson and Ellis 
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[90|. Although the measurement is extremely challeng- 
ing, it may be feasible over a 20 year period by precision 
measurements of the Lyman-a forest in the redshift range 
2 < z < 5 with the next generation of Extremely Large 
Telescopes [9l|. 

For FLRW models one has 

^ = i7„(l + z)-i/(z) (58) 

which in the case of a ACDM model with possible spatial 
curvature leads directly to 

i-^^t = (1 + ^) - V"mo(i + + ^^Ao + ^^.o(i + ^)'- 

(59) 

For the timescape model one has an expression identical 
to (|58p in terms of the dressed Hubble parameter if the 
time derivative is take with respect to wall time, r. Using 
(|B8l) we find that 

1 dz H 

= l + z-^ 60 

H^t {2t + Zbf 

where t is given implicitly in terms of z by (|38p . 




FIG. 12: The function for the timescape model with 

fvO ~ 0.762 (solid line) is compared to H~^^ for three spa- 
tially flat ACDM models with the same values of {^j^jq, ^^ao) 
as in FigM (dashed lines): (i) = (0.249,0.751); 

(ii) (t^^^o'^Ao) = (0.279,0.721); (iii) (J^^.o-^ao) = 
(0.34,0.66). 

In Fig.[T2]we compare H^^^ for the best-fit timescape 
model with /vo = 0.762 to the equivalent function for 
three different spatially flat ACDM models. What is no- 
table is that the curve for the timescape model is con- 
siderably flatter than those of the ACDM models. The 
origin of this feature may be understood qualitatively 
to arise from the fact that the magnitude of the appar- 
ent acceleration is considerably smaller in the timescape 
model, as compared to the magnitude of the acceleration 



in ACDM models. For models in which there is no ap- 
parent acceleration whatsoever, one finds that H~^j^ is 
always negative. If there is cosmic acceleration, real or 
apparent, at late epochs then H~^^ will become posi- 
tive at low redshifts, though at a somewhat larger redshift 
than that at which acceleration is deemed to have begun. 

Fig. [12] demonstrates that a very clear signal of dif- 
ferences in the redshift time drift between the timescape 
model and ACDM models might be determined at low 
redshifts when H~^^ should be positive. In partic- 
ular, the magnitude of H~^^ is considerably smaller 
for the timescape model as compared to ACDM mod- 
els. Observationally, however, it is expected that mea- 
surements will be best determined for sources in the 
Lyman a forest in the range, 2 < z < 5. At such 
redshifts the magnitude of the drift is somewhat more 
pronounced in the case of the ACDM models. For a 
source at z = 4, over a period of St — 10 years we 
would have 6z = —3.3 x 10~^° for the timescape model 
with /vo = 0.762 and = 61.7 km sec"iMpc~\ 
By comparison, for a spatially flat ACDM model with 
= 70.5 km sec^^ Mpc^^ [44] a source at z = 4 would 
over ten years give 6z = —4.7 x 10^^° for (f^jy^Q, ^j^q) — 
(0.249,0.751), and Sz = -7.0 x lO-^° for (f7^^g,rj^g) = 
(0.279,0.721). 

VIII. DISCUSSION 

In conclusion, the combination of tests we have de- 
scribed here have the potential to decide between the 
timescape cosmology, the ACDM cosmology, and other 
homogeneous isotropic cosmologies with other sources of 
dark energy. A number of the tests have been devised 
by other researchers with homogeneous dark energy cos- 
mologies in mind. In these cases, the results of indepen- 
dent analyses performed to date are encouraging for the 
timescape model. In particular, 

• A study of w{z) from recent datasets by Zhao and 
Zhang ^] provides mild evidence at the 95% con- 
fidence level for an effective ^(z) which crosses the 
"phantom divide" near the redshift z ~ 0.46 indi- 
cated in Fig. EJa), with ■w{z) -I- 1 of the same sign 
over the relevant redshift ranges for z < 1; 

• Fits of classes of empirical w(z) functions by 
Shafieloo, Sahni and Starobinsky [80] yield, in 
the best-fit case, an Om{z) function with inter- 
cept Om{0) which appears to coincide with the 
timescape expectation, Oto(O) ~ 0.64; 

• Studies of the BAG scale in SSDS-DR6 data by 
Gaztahaga et al (sT], [s^ yield a relative mass frac- 
tion of baryonic matter to nonbaryonic dark mat- 
ter, which is higher than the WMAP5 expectation 
with a FLRW cosmology, but which is perfectly 
consistent with the timescape model fit to the an- 
gular scale of the sound horizon (38i] . 
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While one can conceive of dark energy models with a w{z) 
which mimics Fig.lJJa) at redshifts z < 1, there is no rea- 
son to expect a different normalization of rt^^/flj^j^ for 
such models. Indeed primordial nucleosynthesis bounds 
are a very strong constraint on all cosmological models. 
It is precisely because the mean CMB temperature at a 
volume-average location unbound to physical structures 
in a void is cooler in the timescape scenario than the 
mean temperature we measure in a galaxy, that a dif- 
ferent normalization of the primordial baryon-to-photon 
ratio relative to present epoch cosmological parameters 
is obtained. This would not be true for any homogeneous 
isotropic cosmology, regardless of the type of dark energy 
fluid. 

Other future tests discussed in this paper also 
have definitive predictions. The expectation for the 
(in)homogeneity test of Clarkson, Bassett and Lu [stI ]. 
yields a diagnostic B{z) which is both distinctively differ- 
ent from FLRW models with spatial curvature as shown 
in Fig. [TOl and from LTB models. The time-drift of cos- 
mological redshifts would be most definitively tested by 
monitoring as many redshifts as possible in the range 
1 < z < 2. As shown in Fig. [121 in this range H^^^ 
should be very close to zero, and only very marginally 
positive as compared to the ACDM expectation. In the 
redshift range, 2 < 2; < 5, which is expected to be the 
range most readily tested with the next generation of ex- 
tremely large telescopes, the function H~-^^ will have a 
flatter z-dependence for the timescape model than com- 
parable ACDM models, as is seen in Fig. [12] The redshift 
range 2 < z < 8 can also be tested by GRB Hubble dia- 
grams, and initial investigations are in progress [TOj . 

This paper has considered tests on scales greater than 
that of statistical homogeneity. There are many other 
such tests in addition to those which we have discussed. 
A number of these involve the CMB, such as the de- 
termination of the amplitude of the late-time integrated 
Sachs- Wolfe effect. Such tests require first a computa- 
tion of the detailed structure of the CMB acoustic peaks, 
recalibrated to the timescape cosmology. This is a very 
complicated task, which is why it has not been attempted 
here. However, it is an important goal for future work. 

Below the scale of statistical homogeneity we expect 
to see apparent variance in the Hubble parameter, with 
a peak value 17% larger than the dressed global average 
value, measured over the scale of the dominant void frac- 
tion of 30h~^ Mpc. Since voids dominate by volume, a 
spherically symmetric average out to a fixed redshift will 
yield generally higher values until we average over vol- 
umes for which a typical line of sight intersects as many 
walls and voids as the global average. That is, the spher- 
ically averaged Hubble parameter should decrease from 
a maximum at the 3Qh~^ Mpc scale to the global aver- 
age value at roughly the 100/i~^ Mpc scale. This general 
pattern is indeed borne out by the analysis of Li and 
Schwarz [s^. Much more detailed predictions of the ex- 
pected variance could be made for the timescape model, 
by performing Monte Carlo simulations assuming a rea- 



sonable distribution of voids and minivoids packed into 
100/i~^ Mpc spheres. This is an important goal for fu- 
ture work, as it would give a Hubble bubble feature with 
unique characteristics, providing a test of a feature for 
which there is no counterpart in the standard cosmology. 

The recent determination of = 74.2 zL 

3.6 km sec-i Mpc"^ by the SHOES survey ^ does pro- 
vide a challenge for the timescape model. However, as we 
have just noted, in the timescape scenario spatial cur- 
vature gradients and apparent variance in the Hubble 
flow below the scale of statistical homogeneity introduce 
systematic issues which complicate the determination of 
Riess et al [93] have done a very careful analysis, 
and make efforts to account for a Hubble bubble - which 
they cut off at z = 0.023, approximately two thirds of 
the scale of statistical homogeneity. However, while they 
do not use supernovae with z < 0.023 in the measure- 
ment of the Hubble flow, their calibration of the distance 
ladder is necessarily made on nearby scales, in particu- 
lar using the maser distance to NGC 4258, at 7.2 ± 0.5 
Mpc, as an anchor. In the timescape scenario the effects 
of spatial inhomogeneity and spatial curvature gradients 
are greatest on scales up to 30h^^ Mpc. Given that our 
own galaxy appears to be in a fllament, this may have 
an impact in calibrating standard candles in the distance 
ladder. 

The megamaser project ^9§\ will therefore provide an 
interesting test, as it will yield purely geometric distances 
- independent of standard candle calibrations - on scales 
much larger than has been tested to date [11]. The rel- 
evant scales are considered to be well into the Hubble 
flow in the standard cosmology, and if distances of order 
Mpc could be measured, would represent a sub- 
stantial up to a large fraction of the scale of statistical 
homogeneity. The expectation in the timescape scenario 
is that provided such sources are sampled in directions in 
which the line of sight passes though a variety of different 
density fields, then there should be variance in the values 
of the Hubble constant so derived. The sample of maser 
distances required to test the statistical expectations of 
the timescape scenario would be considerably larger than 
the ten or so masers currently under investigation, but 
may become feasible in coming decades. 

In comparing future measurements with model predic- 
tions it is important not only to extend the timescape 
model to develop counterparts of all the standard tests 
of the FLRW models, but also to carefully examine the 
methods by which astronomical data is reduced, as in 
many cases the standard cosmology is either explicitly or 
implicitly assumed. As one case in point, BAG analyses 
at present typically use a transformation to Fourier space 
and the use of spectral transfer functions calibrated to the 
FLRW models. Thus while the results of Gaztanaga et al 
[Sll [ssj are suggestive in that they find results in agree- 
ment with our expected H^p/n^p - which is the physical 
parameter responsible for the degree of baryon drag in 
the primordial plasma - in applying results of indepen- 
dent analyses to the timescape model one must exercise 
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caution until each step in the BAO data reduction is un- 
derstood directly from calibrations with the timescape 
model. 

Another important case in which data reduction must 
be carefully considered is that of supernovae. It was 
recently pointed out [qHI that on Bayesian evidence 
the timescape model is disfavoured as compared to the 
ACDM model using the Union [41j and Constitution [i^l 
compilations. However, the Union and Constitution data 
sets have been reduced using the SALT method in which 
one simultaneously marginalizes over both empirical light 
curve parameters and cosmological parameters, assuming 
a FLRW cosmology. Hicken et al [43] discuss and com- 
pare four different methods of data reduction: SALT, 
SALT2, MLCS31 and MLCS17. They find some system- 
atic differences between the methods; for example, the 
SALT methods give larger scatter at higher redshifts. 

As will be discussed in a forthcoming paper [i^ , use of 
the MLCS17 reduced data gives a different picture to the 
conclusions drawn by Kwan, Francis and Lewis [o^. In 
particular, analysis of the MLCS17-reduced 372 Snela 
of Hicken et al gives Bayesian evidence which favours 
the timescape model over the ACDM model. Thus there 
are already enough supernovae in principle to distinguish 
between the models; except that systematic uncertainties 
in the empirical methods by which standard candles are 
standardized at present limit the conclusions that can be 
drawn. Such issues are likely to also be a feature of many 
other astrophysical observations, and thus it is important 
that as many independent tests as possible are devised, 
and carried out carefully in a way in which any model- 
dependent assumptions are scrutinized. 

It is hoped that the tests discussed in this paper 
will provide a basis for comparing the ACDM model 
with a physically well-grounded competing cosmological 
model. To fully compete, much further development of 
the timescape model is of course required. The standard 
cosmology consists of a base model for expansion of the 
universe - the FLRW model dating from the 1920s - on 
top of which a sophisticated superstructure has been built 
over the last few decades. This superstructure includes 
features such as the generation of initial conditions from 
inflation, the bottom-up hierarchical structure formation 
process, and the results of large-scale structure simula- 
tions using Newtonian gravity on top of the base expan- 
sion. The model of refs. [H, [2^ replaces the base expan- 
sion of the FLRW model by an average expansion which 
is not based on the Friedmann equation, and this paper 
has explored a number of tests which can be performed 
based solely on the average geometrical properties. 

Many current cosmological tests of the standard cos- 
mology - including detailed analysis of the CMB, galaxy 
clustering, redshift space distortions and weak lensing - 
can only be extended to the timescape model once the 
standard cosmology superstructure built on top of the 
FLRW model is adapted to the timescape model to un- 
derstand the growth of structure at a more detailed level. 
Although this may seem a daunting task, it is perhaps 



not as quite a tall order as one might at first think. In 
particular, the differences from a standard FLRW model 
with inflationary initial conditions at last scattering are 
negligible, and consequently many large portions of the 
standard cosmology would not change. In particular, 
the mechanisms of physical processes are largely still the 
same, but what does change is the relationship of present 
average cosmological parameters to the initial perturba- 
tions. Rederivation of the standard cosmology super- 
structure may largely be an issue of recalibration. Where 
the calculations involve transfer functions that relate ini- 
tial perturbation spectra to their time evolved distribu- 
tions, such recalibrations may be quite nontrivial, how- 
ever. Thus a careful first principles re-examination is 
required. This is left to future work. 
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APPENDIX A: GENERAL TWO-SCALE 
SOLUTION TO THE BUCHERT EQUATIONS 

The general solution for the two-scale [96| Buchert 
equations pT)) for the independent functions a{t) 

and fv{t) is given implicitly by [29| 
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where u = fv^^^a/a^ — fvi'^^a^/a^^ is proportional to 
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being constants of integration, while /vo, -Hq, 
and f2feo f-r^ the present epoch values of /v, H, Clj^^ 
and rij, respectively. Since /wi^^'^a-w — (l~/v)^^'^a, 
Eq. (jAl|) may also be written as — a^Qt^^^^ where 
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The lapse function, 7, bare matter density parameter, 
fi^^, and void fraction, /v, satisfy the integral constraint 
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Furthermore, H„ = 2/(3t), while Hy = H„/hr where 



Of the six constants e^, t^, /vo, J?q, f^jv/p ^^"^ '^O' 
only four are independent since there are additional con- 
straints [291 
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where the age of the universe in volume-average time is 
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on account of (|Aip . 

Of the four independent parameters, two can be elimi- 
nated by demanding priors at the surface of last scatter- 
ing which are consistent with the evidence of the CMB. 
The redshift of the surface of last scattering relative to 
wall observers at the present epoch, z ~ 1100, is fixed by 
the ratio of our locally measured CMB temperature rela- 
tive to the temperature scale of matter-radiation decou- 
pling and recombination, which is for the most part de- 
termined by the binding energy of hydrogen. We require 
that the velocity perturbations and density perturbations 
at this epoch when z ~ 1100 are consistent with observa- 
tion. For example, we can fix velocity perturbations by 



demanding 1 — hri ~ 10~^ and density perturbations by 
restricting /vi. Physically, /vi is to be understood as the 
fraction of our present horizon volume, TC, which by cos- 
mic variance was in uncompensated underdense pertur- 
bations at last scattering. If this uncompensated fraction 
is viewed as a single density perturbation then 
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We might demand 10 ^,—10 ^}, which means 

we might take /vi G {10^^,10^^}, depending on what 
values of {5p/ p)vi are acceptable for the nonbaryonic dark 
matter power spectrum. 

Once values of hri and /vi are specified, then by (I14p . 
(dropping the index w), = 1 — f^i + fvih~j^, while the 
initial matter density parameter, ^pj^, is fixed in terms 
of 7j, ei and /vi by (jA3[) . At the present epoch, the 
integral constraint (jASp . combined with the relation for 
the cosmological redshift determined by wall observers, 
z+l = a 7/(57 ), gives 
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where we have used (IA4[) to express a^/cii in terms of hri 
and other parameters in the last step. We evaluate both 
(|Aip and (jA2[) at the present epoch and at the time 
of last scattering, ij, and compare them at each epoch to 
eliminate and ti. We then further eliminate from 
the two resulting expressions to also obtain 
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For fixed z^, /vi and hri the combination of eqs. (jA5[) . 
(|A9p and (jAll[) determines three of the parameters 
fjfeOj Eii /vo}i leaving one independent parameter 
in addition to H^. Of course, the values of z^, /vi and 
hri which are consistent with the observed CMB, do 



vary over some small ranges. Given the existence of the 
tracker solution, however, these small variations do not 
significantly affect macroscopic cosmological parameters. 
The macroscopic properties of the universe depend signif- 
icantly on the two independent parameters H„ and fvo- 



20 



APPENDIX B: TRACKER SOLUTION TO THE 
BUCHERT EQUATIONS 

As noted in ref. [2§|, setting = in the gen- 
eral solution gives a solution which is a strong attrac- 
tor in the phase space. Physically this solution repre- 
sents one in which the void regions expand as empty 
Milne universes in volume average time, av= dvot, where 

tion is given by 



Q'''n H r^'^f .-1/3 



and hr — 2/3. The solu- 
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with two independent parameters and /vq. 

All other quantities of interest may be determined from 
(IBT|) and (|52|l . For example, the parameters P^ - ipH]) 
are given by 
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where 6 = (1 - /vo)(2 -h /vo)/[9/voi?o]' in ([39]). From 
(|B3|) -(|B5 P we obtain equivalent expressions for their 



present epoch values, ^j^q, ^ko and ^Qf^, in terms of /, 
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or (2 + /vo)/(3iJp). The bare Hubble parameter, 

lapse function and dressed Hubble are given respectively 

by 
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For a number of the tests described in the paper, it is 
necessary to perform derivative of observational quanti- 
ties with the respect to the redshift, z, as given by (|38|1 . 
For this purpose, a useful intermediate step is provided 

by 
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which follows from (j38|l . Thus, for example, by ((38 
and (|BT0)) . 
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The above expressions all involve volume-average 
time, *. To relate them to wall time, r, which is as- 
sumed to be a good approximation to the time measured 
by typical observers in galaxies, one has to invert the 
relation 
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where fi^^g = 5(1 ^ /vo)(2 + /vo) is the present epoch 
dressed matter density. 
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